“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1975 


An alternative to nonlinear estimation in gun 
fire control systems 


Brands, Robert Keith 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/20721 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


, 






/ - | 
7 > —_ 


AN ALTERNATIVE TO NONLINEAR ESTIMATION 
| IN GUN FIRE CONTROL SYSTEMS 





| 


| | 
| 


‘ 





| 
| 
ee | Robert Keith Brands 


= = “a — = 





we 


‘. 








» AN ALTERNATIVE TO NONLINEAR ESTIMATION 
IN 
GUN FIRE CONTROL SYSTEMS 


by 
Robert Keith Brands 


December 1975 





Sl ie eles Gag SOT. Sime a 


Approved for public release; distribution unlimited. 


T17166¢ 





Wimelkasc it led 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE BEFORE COMPLETING FORM 


REPORT NUMBER 3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subdtitie) S. TYPE OF REPORT & PERIOD COVERED 


t ° 
An Alternative to Nonlinear Estimation Master's Thesis ; 


: ; December 1975 
Mmcui fare Control Systems 
6. PERFORMING ORG. REPORT NUMBER 









READ INSTRUCTIONS 


2. GOVT ACCESSION NO. 








7. AUTHOR(a@) i 8. CONTRACT OR GRANT NUMBER(A) 
ein nas 
CEE rand 
SG / 2 : IN NAME ANO AOORESS 10. PROGRAM ELEMENT, PROJECT, TASK 
Ee ons ak, AREA & WORK UNIT NUMBERS 
jate School 
"Ornada 95940 
JLMo Bye G+ dl ! 
ME ANO AOORESS 12. REPORT OATE 
= ; December 1975 
Ge Pen 0% c late School 13. NUMBER OF PAGES 
‘fornia 93940 
a 4 enn 3) ae 
E& ADORESS//f different from Controlling Office) 15. SECURITY CLASS, (of thie report) 
late School Unclassified 
Se Z Ye, ‘ornia 93940 
; , 1Sa. DECLASSIFICATION/ DOWNGRADING 
SCHEOULE 


(of thie Repcrt) i 


DpUMICe release; distribution unlimited. 


(of the ebetreact entered {fn Block 20, if different from Report) 


18. KEY WORDS (Continue on raverce eide if noceceary and identify by block number) 


Kalman Filter bineareestimator 


Spherical Coordinate 


Fire Control System 





20. AGSTRACT (Continue on reverae efde if neceersary and identify by block number) 

A Kalman filter with assumed linear system dynamics was 
derived from a spherical coordinate state vector and applied 
to a fire control system tracking environment. Two models 
of target dynamics were developed: a constant-velocity 
Nedelandeamacerrelated random acceleration model. The 
predicted position performance of the filters derived from 
the two basic models was optimized using Monte Carlo methods 


DD . en 53 1473 EDITION OF 1 Nov 6815 CASOLETE Wines ics od 
(Page 1) S/N 0102-014- 6601 | 


etter eg wg a te PP PC 
SECURITY CLASSIFICATION OF THIS PAGE (Mhan Data Entered) 





Wrelassitied / 
SECURITY CLASSIFICATION OF THIS PAGE(When Deta Entered) 


meowtistsa set Of teSt trajectories that include constant- 
velocity and maneuvering target profiles. The results of 
the Monte Carlo simulation of these filters are compared 
Wiha tne results Obtained from a filter derived from a 
Cartesian Coordinate state vector. Switch-on-range 
adaptive filters were developed from the two basic models 
and evaluated by Monte Carlo methods. The results of 

the simulations of the adaptive filters with the two dif- 
erent State vectors are compared. 


DD For 1473 ACK) ie 
m (B Unclassified 


an 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 








Am ALterpnative teoemNenlinear Estimation 
iat 


Gun Fire Control Systems 


by 


Robert Keith Brands 
Lieutenant, United States Navy 
Beow,eornedd State Collece, 1969 


SUbDNTtued in Dartaal tultvililment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 





ABSTRACT 


A Kalman filter with assumed linear system dynamics was 
derived from a spherical coordinate state vector and applied 
moma tire Control system tracking environment. Two models 
of target dynamics were developed: a constant-velocity 
model and a correlated random acceleration model. The 
predicted position performance of the filters derived from 
the two basic models was optimized using Monte Carlo methods 
against a set of test trajectories that include constant- 
Mmoctty Maneuvering target profiles. The results of the 
Pemee) Carlo simulation of these filters are compared with the 
results obtained from a filter derived from a Cartesian 
Meenainate State vector. Switch-on-range adaptive filters 
were developed from the two basic models and evaluated by 
Monte Carlo methods. The results of the simulations of the 
adaptive filters with the two different state vectors are 


compared. 
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I. INTRODUCTION 

A. AN ALTERNATIVE TO NONLINEAR ESTIMATION IN GUN FIRE 

CONTROL SYSTEMS 

The increased performance characteristics of modern 
aircraft and missiles has made apparent the need for more 
effective fire control systems. With the introduction of 
digital fire control computers, Kaiman filtering tech- 
niques offer an attractive means of improving the effective- 
ness of a fire control system. Indeed, several recent 
Studies conducted at the Naval Postgraduate School [Refs. 
1-6] have specifically addressed the application of Kalman 
ieecrincg techniques to fire control systems. It is 
interesting to note that the references cited above all have 
one distinctive feature in common; they use Cartesian 
Meencinate state vectors. A Cartesian Coordinate state 
Memror offers the attraction of linear plant dynamics for a 
memouant-velocity target but has the disadvantage of non- 
linear measurement equations when used in a fire control 
System that measures the target's position in spherical 
Semrainates. One possible solution to the nonlinear measure- 
ment problem is to process the sensor measurements through 
a pretilter as in Figure 1. The output of the sensor 
consists of noise corrupted measurements of the target's 


Mesation in spherical coordinates 


iy 
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S, (K) 6 (alee n, (CK) 


s,(K) = ¢(K) + n,(K) (1.1) 


Sz (K) nek) + nz (K) 


which serve as inputs to the prefilter. The appropriate 


coordinate transformation equations from Figure 2 are 


Cie = rik) sind ( Keosd (kK) 
y(K) = r(K)cos¢ (K)cos6 (Kk) (lez) 
z(K) = r(K)siné(K) 


where it has been shown [Ref. 3] that v(K), the measurement 
noise process presented to the tracking filter is a non- 
linear, position-dependent function of the variances of 
Mmiemrange, bearing, and elevation errors. Clearly then, the 
measurements Wecone correlated in the coordinate conversion 
meeeess. Full implementation of such a filtering scheme 
could conceivably involve much on-line computation since 
the statistics for v (K) are no longer stationary. Another 
fereinct disadvantage of the prefiltering method lies in 
the errors generated due to the added finite precision 
Moichmetic of the coordinate conversion process. Such 
errors may be qualitatively described as a noise process 
Pno-c Statistical properties are difficult at best to 
identify. An excellent treatment of this specific topic is 


found in Gelb [Ref. 7]. 
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LMOtleredpprodacn to the nonlinear measurement process 
is to linearize the measurement equations about some 
reference or nominal trajectory. Unfortunately, such an 
approach is usually accomplished by the implementation of 
emeextended or iterative extended Kalman filter with the 
merendant increase in On-line computing requirements. 

Consider the filtering problem with a spherical 


coordinate state vector 
_ e te e ee e ee T 
Oana Ere 65.0, 0 | (1.4) 


Such a state vector has the advantage of linear measurement 
equations eliminating the correlation and position dependence 
of the measurement noise process. The disadvantage of a 
Smeraical State vector lies in the coupled and nonlinear 
System-target-dynamics, even for a constant-velocity target. 
Once ein, the implementation of this filter would normally 
be with an extended or iterative extended Kalman filter. 

mem tnesis proposes a filter having a spherical coordinate 
State vector with assumed linear system dynamics. In view 
of the fact that any system model is an approximation to a 
physical system, this seems to be a reasonable assumption. 
Implementation of this would be a great deal simpler than the 
previously mentioned filters. A filter designed on this 
basis would then of necessity be sub-optimal. However, 

given acceptable performance, the sub-optimal filter might 


Pespreterred to the optimal nonlinear filter. Faced then, 
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with a sub-optimal filter due to modeling inaccuracies, 

a technique must be found to compensate for the unmodeled 
nonlinear properties of the system. Jazwinski [Ref. 8] has 
proposed two methods for the compensation of inaccurate 
mepriori Knowledge of the system. The first involves a 
fictitious noise input to the system, while the second 
Overweights the most recent data. The net result of either 
method is an increase in the filter bandwidth, making the 
molec, more responsive to system transients but increasing 
its sensitivity to the measurement noise process. The 
fictitious noise input method was chosen as the most promis- 
ing and was used in this thesis. So that some quantifiable 
Heestre of the feasibility of these filters is possible, 
this thesis paralleled Parr's work [Ref. 4] and comparisons 


aves made with his results. 


B. LINEAR FILTERING THEORY 
Mies system Studied in this thesis 1s composed of a 
target and a sensor which gives noise-corrupted measure- 
Memts OL the target's position in spherical coordinates. 
The system may be approximated by the linear difference 
equations 
x(K+1) = $(K+1,K)x(K) + P(K)w(K) 


felieee5:) 
ZAC te) eek) = vi (K) 


The statistics assumed for the stochastic processes w(K) 
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and v (K) are 


E[w(K)] = E[v(K)] = 0 for 


Elv(k)y (J) = R(K) 6, 


Efw(k)w (5) = Q(K)S,, 


E[w(K)v’ (K)] = E[v(K)w’ (K)] 


E[w(K)x'(0)] = E[x(0)w' (K)] 


E[v(K)x(0)] = E[x(0)v  (K)] 


all K Gee) 
Cle) 
GU 
EO ecoy all unk (1.9) 


H 
a 


roe BILL el 0) 


for all K (1.11) 


I 
a) 


Miemlinear, discrete Kalman filter is characterized by the 


difference equations 


X(K/K) = X(K/K-1) + G(K)[z(K) - H(K)x(K/K-1)] qi) 


WMeeirtne One-Step prediction given by 


x (K/K-1) = ¢(K,K-1)x(K-1/K-1) 


(aes 


The Kalman Gain matrix, G(K), is found by solving the 


following equations 


G(K) = P(K/K-1)H? (K) [H(K)P(K/K-1)H (K) + R(K)]” 


P(K/K) = [I - G(K)H(K)}P(K/K-1) 


nel) 


Gai 


P(K+1/K) = $(K+1,K)P(K/K)¢" (K+1,K) + 


POE CS 


(1.16) 


To start the gain calculations, the matrix P(K/K-1) must 
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be initialized. Taking the standard approach, P(0/-1) 


is given by 


T 


P(0/-1) = E[(x(0) - X,)(x(0) - X))°] (1.17) 
where 
E[x(0)] = X, (1.18) 


mirenblased eStimator 1S obtained if the filter is started 


Dyeeset ting 


x(0/-1) = Xp (1.19) 
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tev eeoruEnd OF THE FILTER MODELS 


A. CONSTANT-VELOCITY (1/82) MODEL 

The first filter tested was based on a dynamic system 
model that assumes linear constant-velocity target motion 
imeach Of three coordinate directions. 

The state transition matrix > (K+1,K) and the random 


foeetne Input distribution matrix I'(K) are 


@(K+1,K) = (2.1) 


and 
y 
T 
x 0 0 
T 0 0 
2 
T 
0 x 0 
Tr (K) - (2257) 
a 0 T 0 
2 
T 
0 0 7) 
0 0 T 


Where T is the sampling period. 
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SiMcCmiENeCmottLemVcCeLOr 15 expressed in spherical 


Beondinagtes, the measurement process is linear and is given 


by 
z(K) = H(K)x(K) + v(XK) (225) 


The measurement noise processes for the individual 
Goordinate directions are assumed to be independent 
meoeeesses. The covariance of measurement noise matrix 


Imi) is diagonal and is given by 


rom SS) 
6 (Se) 


R(K) = O (24>) 


Gore 


Because of its simplicity the constant-velocity model 
cannot realistically be expected to perform well against 
all types of target trajectories and may not perform accep- 
tably against maneuvering targets. The performance of the 
constant-velocity model filter can be improved by the use of 
the system random input covariance matrix Q(K). Indeed, 
eeeeot the functions of the Q(K) Naeix Se tO —cCONpensia te 
mon the inadequacies of a simplified model, and for this 
meeltic thesis, Q(K) CHismesoOuparerailily account tor the 


System nonlinearities not explicitly included in the 


system model. The subject of the sensitivity of the 


Zi 








meter '’s N-step prediction performance to the Q(K) matrix 
forms the basis for the simulations of Chapter III. 
Initialization of the filter requires a priori knowledge 
Geetne Statistics of the system initial conditions so that 
@iemegain calculations may be started. For any fixed coordi- 
nate system, it is reasonable to consider the mean of all 


possible system initial conditions as zero or 


B[x(0)] = X, = 0 (2.5) 


ieeras been Shown [Refs. 7 and 9] that under proper 
tmimetalization, the filter will remain unbiased for 


me, 1.e., 
E[e (K/K)] = 10 Ore BIL Is (20) 
where 
e(K/K) = x(K/K) - x(K) Ca?) 
iiesunbiased condition is satisfied if 


x(0/-1) = (2.8) 


Pes 


Initialization of the covariance of estimation error 
matrix requires a somewhat different approach. Following 
Parr [Ref. 4] 1t is quantitatively shown that if the 
diagonal elements of the initial covariance of estimation 
error matrix are made very large, the net effect is that 


the first position estimate is made identically equal to 


a2 








the first system measurement. In this manner the system 
iemeorced into tracking the target rapidly. For all simu- 


lations conducted then, it was assumed that 


Faia) = 8) 
10 


a ae 


L 


B. CORRELATED RANDOM ACCELERATION MODEL 

The constant-velocity model discussed in the previous 
section views any target maneuvers as small perturbations 
about a nominal, constant-velocity trajectory. This tacitly 
assumes that any target maneuvers are uncorrelated in time, 
or from one measurement to the next. Target maneuvers will 
generally take place over a period of n measurements and the 
iomeuae Of mn is dependent on the level of the target's 
Maneuver, 1.e., for a slow turn n is relatively large 
whereas for an evasive type maneuver n is smaller. The 
fieemtnat the accelerations occur over a finite period of 
time correlates the target's accelerations (maneuvers) in 
mies One approach to modeling the effects of time- 
correlated target maneuvers is to use a constant-velocity 
plant driven by correlated noise obtained from a low-pass 


meter of the form 


1s) = (s * a) 20) 
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driven by white noise. The constant a is interpreted as 
the time constant or degree of time correlation of the 
target's maneuver. 

For one dimension, the overall model transfer function 
1s 
te) 


T(s) = —-———-—— 
se + 5 C ) 


(2eias 
femacpicted in Figure 3. 
bor the three-dimensional case, the filter is ninth 


@raer and has the state vector 
» ® ee e ee e ee dp 
x(K) = [r,7,7,¢,6,6,0,0,6] Car 


The correlated random acceleration model is discussed in 
detail by McKinley [Ref. 5] and Singer [Ref. 10]. From the 


m@ereu references, the (3 x 3) state transition matrix 


2505 and the (3 x i) random forcing input vector lr, (T) 


are for the single dimension with assumed linear system 


dynamics 
fk it -~aT 
a0) E Aon —= (1 ae ) 
co 
6.(T) = | 0.0 1.0 See?) Ss) 
0.0 0.0 Qo | 
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and 
ae 
2 a o 2 
T.(T) = to 2 ee Cai 
~1 0; 
(lame 5) 


For the ninth-order filter the state transition matrix 


o(K+1,K) is formed by 


g,(T) 0 L 
eck 1K) = 0 dp (T) 0 (215) 
" alls 
where o,(T), dp(T), and op (T) are individually given by 
>; (T) in Equation (2.13). In a similar manner, for the 
ninth-order filter the three-dimensional system matrix 
r(K) is given by 
0 0 
r¢(K) = P, (T) 0 (2G) 
0 rg (Ty 





where al; P,(T),° P.C) are individually given by a) 


in Equation (2.14). 
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The covariance of measurement noise matrix R(K) 1s the 
same as given in the preceding section by Equation (2.4). 
feetalization of the filter is carried out in exactly the 
same fashion as described in the previous section by 
Memations (2.5), (2.8), and (2.9). As in the preceding 
section, Q(K), the covariance of the random forcing input 
white noise process presented to the coloring filter, was 
used to compensate for the unmodeled dynamic nonlinearities. 
iMaemsensitivity of the filter's N-step prediction perfor- 
mance to the Q(K) matrix selected is treated in 
Chapter III. 

In the interest of Saving computer time, the correlated 
random acceleration model was tested using a single value 
for a, the maneuver time constant. Since Parr [Ref. 4] 
reported that the best overall performance in his simula- 
tions was obtained from a value of a = 0.5, that value was 
used in the study of this model. Additionally, the maneuver 
time constant a, was fixed at the same value for each of 
the three coordinate directions, though this constraint 


need not necessarily be imposed on the system. 
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Bite om owOm fie MONTE CARLO’ SIMULATION 
CPi ROLOSkE be LINEAR FILTERS 


Poe EXPERIMENTAL PROCEDURE 

Monte Carlo simulation techniques were used in the 
optimization and testing of the filters developed in 
Chapter II. A listing and description of the FORTRAN 
program used is given in Appendix A. All results are the 
products of 100-member ensemble simulations. Test trajec- 
tories 1-8 were used in the optimization and testing process 
moeeepocndaix B gives a description of the profile for each 
MmemectOry. the profiles of trajectories 1, 2, 3, and 4 are 
identical to the short-range constant-velocity, short-range 
Womemvering, long-range constant-velocity, and long-range 
maneuvering profiles, respectively, used by Parr [Ref. 4]. 
Trajectories 5-8 were included to examine the filter's 
performance against more Severe target maneuvers. 

The elements of R(K), the covariance of measurement 


noise process matrix, were fixed for all simulations at 


Op = Seu vas: 
Op = 0.01722 Deg. (G3z71)) 
o, = 0.01722 Deg. 


fie values given in Equation (3.1) are those used by Parr 


[Ref. 4] and are used in this thesis in an attempt to 
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duplicate Parr's testing environment, and minimize procedural 
(moerepancies. 

The performance evaluation of the filters developed in 
Chapter II is based on their N-step prediction or position- 
Sate prediction accuracy. A detailed description of the 
method used in the determination of the error in the 
Predicted position is given in Appendix C. The results 
Peesented in this chapter are the percentages of “hits" 
ana, ''scares' comptited oWer the entire ensemble. A “hit' 
is defined as an error in the predicted position of less 
Giaeetont yvards, while a "scare'' is Wefined’as am error 
in the predicted position of less than 16 yards, but ¢greater 
than eight yards. 

Initialization of each filter simulated was identical 
Mmeati® trials and Ys described’'in Chapter [I by Equ@tions 
wey, (2-8), and (2.9). 

The elements of the covariance of the random forcing 
input matrix Q(K) were determined by a two-dimensional 
optimization search. The assumption of linear uncoupled 
dynamics eliminates any spatial correlation, resulting in 


a diagonal Q(K) matrix of the forn 


q,. 0 0 
Quk) = [0 ls 0 (3.2) 
10 0 qe | 
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If the system dynamics were truly linear, then the diagonal 
elements of (3.2) would represent the variances of the 

random forcing input to the system. However, since Q(K) 

is also used to compensate for inaccuracies in the system 
model, dy> Wp and dp must be adjusted accordingly. To 
faerity the optimization process for both models, dp and dp 
were constrained to be equal throughout. The two-dimensional 
Search optimized the filter's poSition-state prediction 


accuracy as a function of two variables, q.. and the pair 


ii 
Spee tp: 
In the case of the constant-velocity model, the elements 
of the matrix Q(K) serve as a means by which the bandwidth 
of the filter is increased to ensure that the filter will 
Tespond to target accelerations. This, of course, is in 
addition to compensation for the unmodeled nonlinearities of 
the system dynamics. Because of the dual role played by the 
matrix Q(K), re force search was done to determine the 
elements of Q(K) that gave the best position-state prediction 
performance. In a similar manner, for the filter based on 
the correlated random acceleration model, the bandwidth 
must be increased to permit sufficient response to targets 
exhibiting non-zero acceleration rates. Here again brute 
force search techniques were used to determine the Q(K) 
matrix that gave the best position prediction performance. 
For both models, once the optimum Q(K) matrix was 


i@entitived £o60r a specific trajectory, that particular 
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filter configuration was tested against the remaining 
trajectories. The net result being that each optimized 
filter was tested against seven additional trajectories to 
ensure that the optimum Q(K) had been found, and to deter- 
mame the effect of using a Q(K) selected for a track of 
one type against targets having different maneuvering 


@meracteristics. 


B. SIMULATION RESULTS FOR THE CONSTANT-VELOCITY (L/S“) MODEL 
ie results presented in this section are for the opti- 

iweeecad tilter configurations. The term "Filter 1'’ denotes 

the filter with the specific Q(K) matrix that gave the best 

Peetticon prediction performance against trajectory 1. Thus, 

miemsimulation results for eight filters, or more precisely, 

eight filter configurations are presented. 

Peepointed OUt in the preceding section, the assumption 
Gmerinear uncoupled dynamics results in the diagonal Q(K) 
matrix given by Equation (3.2). Table I gives the diagonal 
elements of the Q(K) matrix that provided the best position 
Peee@reti1Oon periormance for each test trajectory. The rela- 
tive magnitude of the diagonal elements for the optimal 
filters are as expected. That is to say, the elements of Q(K) 
iO the non-maneuvering trajectories (1, 3, 7, 8) are lower 
in magnitude than those for the maneuvering trajectories 
(2, 4, 5, 6). Similarly, the elements of Q(K) for the long- 
range maneuvering trajectories are lower in magnitude than 


those for the short-range maneuvering trajectories. Similar 


SI 








Table I 


DIAGONAL ELEMENTS OF THE OPTIMAL Q(K) MATRICES 


FOR THE CONSTANT-VELOCITY (1/S%) MODEL 


FILTER 011 - 
1 1.00 x 107° 8.58 x 10°” 
: De Os Oleleos 
3 1.00 x 10°° eS x 19% 
4 0.025 0.006553 
5 0.0075 0.1372 
6 1.00 x 5.85 2 a | 
/ 0.01 $50 10! 
8 0.005 recut ad 
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.66 x 10° 
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observations may be made for the moderate level maneuvers 
as opposed to the evasive level maneuvers. Note that the 
elements of Q(K) for the long-range evasive maneuver profile 
(trajectory 6) are smaller than those for the long-range 
moderate level maneuver profile (trajectory 4). This appears 
to be a reasonable phenomenon in light of the fact that 
overweighting the estimates for a faster more highly maneu- 
vering target at long range could result in larger errors 
in the predicted position. It follows then that the band- 
width of the optimal filter for the long-range evasive 
maneuver should be somewhat smaller than the bandwidth for 
the moderate level maneuver at long range. 

fable Ji gives the simulation results for the optimal 
fivcer configurations as well as for each optimal filter 
SomeneuratioOn against the othér trajectories. The results 
mieowm are not Surprising in that for a specific trajectory, 
the optimal filter configuration gave the best results. In 
a few isolated cases, other filter configurations gave 
results that approach those for the optimal configuration. 
Mivomen this was not expected, the differences are not con- 
sidered small enough to justify any additional effort in 
the optimization process. 

Table III compares the hit/scare percentages from the 
Simulations of the two models using the two different state 
Mectors. The results for the spherical coordinate state 


vector represent the optimal performance against a specific 


33 





0°0 


0°0 


0°0 


ev) 


56.0 


Ca tees 6 


98 0 


OT 0 


SIBIS2 


Hy ASO OL ieae I 


0°0 


0°0 


0°0 


0° 0 


eu) 


0°0 


i 


0°0 


1TYy 


CG GE 


ou Ce 


Cn ac 


Se 


ie ok 


Vie 


Gs) 


Oa CC 


IIBISe 


8o° OP 


vc OP 


a 0 


0°0 


OTS 


[80 


eae) 


Tc OP 


1TUY 


¢ Ato VOel eal 


Sud Ce adh eee Fae INV ESNOo 


Se ALTE 


oy eb 


ee 


99° RC 


ASO 


sy 1! 


SPY 28 


sie 1h ih 


SIBIS2 


2 Azozoolery 


€0°6 


CU 6 


CU 6 


78° vd 


rosea O¢ 


CU 6 


GS Se 


C02 6 


1TUY 


BO. Vv 


Loy 


89°90 


Omen 


v8°S 


T6'? 


8b OT 


Oilers 


SIeBISe 


€9°98 


69°98 


(Gee sie: 


Vee 


C3 s6 


98°98 


Bite yo 


68°98 


1TYy 


T Ato ..eleAL 


GaZINTidO dai uOd SLIASAY NOIIVIANIsS 


II 91921 


34 


rol Si Ohi Lal 





Pore! 


OZ TE! 


EG it 


ec. 


Se Le 


25 


Rel 


9¢° TT 


ILBISY 


02°69 


nO 69 


cb°89 


SiGe V 


4 BIST 


60°69 


EON Ree 


Oh IL tee) 


ITUS 


g ALOZIOLECIT 


jig eae 


Oars 


92°S 


CL, UC 


€9°8 


SES 


Hee bere 


Sos 


ILEIS 4 


vt'98 


cc 98 


piss 


ood 


0a Cs 


COIS 


08° 19 


GG. 8 


ITUS 


Pe acoavaye Rasp 


0°0 


sw 


0°0 


8) 


pt oO 


fips se | 


wih 


ILBIS% 


oO) neat aie) ene if 


00°42 


96°9 


ct 8 


b0°0 


00°2 


ov 0 


Oat 


4TUS 


Ca oak 


poe Oil 


CO OT 


1a er 


Oc 


OZ. Sat 


vo TT 


Ss 


DLVIS2 


Soe 


20. “VE 


ie ee | 


Fo Lae 


<V[eymateré 


SV CE 


CO VLC 


Ge Cee 


+TU% 


Gg AloqZoeleIy 


(P,3U0D) SUALTIA TSQOW ALIDOTIA- LNVLSNOD 


CiZzIWildoO ant uCias Lins dd NOLL INS 


Il 9TdeL 


volisilk Let 


Ge) 





Gab 2 fell Ore eb BOC a! 


a0 ig 1 SC ae ac 


SLledS% 3TY% s1ed9S% YTYE 


FONVU-ONOT a NV = aero 


SHTYOLOYCVaL ONTYFANANVN 


Lane Lv 00°¢ v°88 
pI°t2 T8°OVv 00°S 68°98 
OLBOS% 2TYs OLBOS% ATYS 
TONVY-ONOT dqUNVd-LdOHs 


Sod Olan Vale whe Taye LN LSNOO 


SCO al OOtdee lt yloNOos dO) aby ts 
TLVNIGUOOD NVISALUVO GNV TVOINAHdS AHL YOa SLTNSAY AO NOSTUVdWOD 


Piles her 


ALVNIGYOO)D 


NVISHLYVO 


TLVNIG@YOOO 
WVyOdands 


YOLIAA ALVIS 


36 





meajyectory. The results tabulated for the Cartesian 
Boeordinate state vector are the best of the results reported 
byerarr [Ref. 4]. 

Tables IV and V show the tracking performance of each 
Optimal filter against its specific trajectory. For 
example, the tracking performance of filter 1 against tra- 
jectory 1 is given, and so on. The tabulated data represent 
Maeemcans of E€Stimation error, defined by Equation (2.7), time 
averaged over the entire trajectory, and the RMS value of the 
deviation of the means of estimation error about their time 
average. 
Geel MULATION RESULTS FOR THE CORRELATED RANDOM 

ACCELERATION MODEL 

The results presented in this section are for the optimized 
mieeer COnfigurations. Since the optimization search for the 
ninth-order filter required a great deal more computer time, 
this model was only optimized against trajectories 1-3. As 
mmc preceding section, the term "Filter 1A" denotes the 
filter with the specific Q(K) matrix that gave the best 
position prediction performance against trajectory l. 

The tabulated data of Table VI gives the elements of 
the Q(K) ldtaces Enateprovided: the soOptimun perrormance 
maaenst each test trajectory. The data appear to be con- 
sistent within the tracking environments presented to the 
merter. Notice that in each configuration the elements of 


the optimum Q(K) matrices for the correlated random 
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Table VI 


DIAGONAL ELEMENTS OF THE OPTIMIZED 





Q(K) MATRICES FOR THE CORRELATED 
RANDOM ACCELERATION MODEL IUD I ees 


Bee R 


471 422 433 
1A 1.00x107” 5.00x10°9 5.00x10 9 
2A 0.100 14.599 14.599 
3A  1.00x10°° 1.46x1077 1.46x10~7 
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acceleration model are larger in magnitude than those for 

the corresponding entries of the constant-velocity model. 
This seems reasonable because this filter is required to 
respond to a transient (due to the non-zero acceleration 
rate) that is expected to be smaller in magnitude than the 
transient encountered by the constant-velocity model filter. 
For this reason the bandwidth of this filter must be in- 
creased to ensure sufficient response to the non-zero 
acceleration rates. During the optimization search, it was 
observed that, with one exception, this filter model is 
extremely sensitive to variations of elements oy) and dz3 
corresponding to acceleration rates normal to the Line-of- 
Sent. this is to be expected, since such accelerations 
cause the greatest displacement from the nominal trajectory 
adding a source for potentially large errors in the pre- 
@iemea position. The single exception noted was for tra- 
Teerory |}. In this specific case the filter was surprisingly 
insensitive to variations in the elements of the Q(K) Met doce 
and its predicted position response was fairly flat over a 
four-decade range Gig WOO. The elements reported 

for Filter 1A are therefore somewhat arbitrary, lying in 

the center of the flat region. Values for the pair G52 = 133 
above and below this région produced inferior results. A 
feeined search grid in the flat region failed to produce 


moe, SlLoniticant change in the results reported for 


Rilter JA. 
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Table VII gives the simulation results for each optimal 
filter configuration and for each optimal filter configuration 
against the two remaining trajectories. Table VIII compares 
the hit/scare percentages from the simulations of the models 
using the two different state vectors. The data presented 
mepeene Spherical coordinate State vector represent the 
Pomemnum performance for the trajectories considered. The 
results reported for the Cartesian Coordinate state vector 
model represent the best of the results reported by Parr 
[Ref. 4]. Tables IX through XI show the tracking per- 
formance for each optimal filter configuration against its 
eeetric trajectory. The results shown are the time- 
averaged means of estimation error and the RMS value of the 
deviation of the means of estimation error from their 


time average. 
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POSITION ESTIMATION PERFORMANCE OF THE 


CORRELATED RANDOM ACCELERATION MODEL 


PARAMETER 
e.. (yds) 


Geeeecy ds) 
Sr 


ey (deg) 
he 


e, (deg) 


er, (deg) 


1A 


87x10. 


pa Vs 


aA Gaels) 


~55x10> 


RSAC 


~43x10> 


Table IX 
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FILTER 


2A 


6.59x10. 


D Se 


2a Selon 


1.40x10. 
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Table X 


RATE ESTIMATION PERFORMANCE OF THE OPTIMIZED 


CORRELATED RANDOM ACCELERATION MODEL FILTERS 


PARAMETER 
eancy/ S) 
%. Wy/s) 
e; (ds) 
< (d/s) 
e4 (d/s) 


Craea/s ) 
fo 2 


FILTER 

1A 2A 

3.77 3.83 

On 29.70 
2.20x1074 1.64x10° 2 

1.76x107° 0.109 

-1.58x107* = -1,. 50x1974 
3.31x10-° 8.14x107° 
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Table XI 


mecuboRAlLON ESTIMATION PERFORMANCE OF THE OPTIMIZED 


CORRELATED RANDOM ACCELERATION MODEL FILTERS 


PARAMETER 
ex (y/s*) 
ag. (y/8") 
es (d/s*) 
o, (a/s") 
eg (d/s*) 


2 
o-,(d/s”) 
“9 


tebe k 
lA IS 3A 
-0.546 -0.549 -0.526 
Bo o2 5259 S17 
ees) i Peal “5 ean = 
as me 
DO Ga kD Gees ik 1759x110 
-4 -4 
-9,70x10 0.49 5. 359x10 
az a5 
Heese eG 0.10 PaO eG 


47 





IV. ADAPTIVE FILTERING 
ee EERIE NG 


A. THE ADAPTIVE BANDWIDTH FILTER 

A single fixed-parameter filter configuration cannot 
Bealistically be expected to perform wel] again Sieh] 
trajectory profiles likely to be encountered. One way of 
approaching this problem is to use a filter that adapts 
its characteristics to 4 Changing CHACKisae environment, 
©.8., a target maneuver is detected. The design problem 
now becomes two-fold; under what circumstances, or when, 
Should the filter's Characteristics be changed and by how 
much should the filter be Changed? Price and Brown [Ref£. 
11] discuss some interesting approaches to the design of 
adaptive filters; unfortunately their complexity is beyond 
the defined Scope of this thesis. 

The single method to be cOnSidered here will be qa 
Simple SWitch-on-range scheme. The motivation for such a 
Scheme lies in the Simplicity of its implementation. 
Additionally, Parr reported [Ref. 4] that of the methods 
he investigated, the Switch-on-range method provided the 
best overall performance in his Simulations. When filter 
iercking is initiated, a low Q(K) gain schedule (corre- 
Sponding to a narrow bandwidth filter) is used to provide 
the optimum performance against long-range targets. When 
the target's measured range reaches a predetermined 


SWitching point, the filter switches to a high Q(K) gain 


48 





schedule (corresponding to a wideband filter) for the 
short-range tracking situation. 

Both filter models developed in Chapter II are simulated 
in a switch-on-range adaptive mode of operation. Again, 
to duplicate Parr'’s testing environment, a switching range 
of 3333.33 yards is used, and the filters are tested 
Meet trayectories 9 - 14. The profiles for trajectories 
tmeeere are identical to Parr's trajectories 1 - 6. The 
results reported for the adaptive simulations are the pro- 
ducts of 80-member ensemble Monte Carlo simulations. For 
micesamulations of the constant-velocity model, the low 
Q(K) Painiesenecanle ised was that irom Filter 3. The data 
fmeapte Il show that Filter 3 gave the best overall per- 
bemmance against the long-range trajectories treated in 
Chapter III]. Similar considerations yield Filter 4 as 
the source for the high Q(K) gain schedule for the short- 
range tracking environment. When filter tracking is 
initiated, the Kalman gains used are those for Filter 3. 
When the target's measured range becomes less than or equal 
to the switching range of 3333.33 yards, the Kalman gain 
schedule is reset to the steady-state gain values obtained 
from the gain schedule for Filter 4. Similar consideration 
of the data presented in Table VII gives Filter 2A as the 
source Of the high Q(K) Paineseneauhe swhtske Filters 5A 
provides the low Q(K) gain schedule. As in the above 


simulation, the low Q(K) gain schedule is used until the 
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mmeeenine range (5553.35 yards) is reached. At that point, 
the Kalman gain schedule is reset to the steady-state 
Meaimevalues from gain schedule for Filter 2A. Initializa- 
tion of the filters was the same for all simulations con- 
ducted and is described by Equations (2.5), (2.8), and 
(2.9). The statistics assumed in all simulations for the 
measurement noise process v (K) are given by Equations (1.6) 
amide (3.1). 
B. SIMULATION RESULTS FOR THE CONSTANT-VELOCITY MODEL 

OWITCH-ON-RANGE ADAPTIVE FILTER 

The constant-velocity model filter was tested ina 
SWitch-on-range adaptive mode of operation against trajec- 
tories 9 - 14. Table XII shows the filter's position pre- 
dietion accuracy results. The data tabulated in Table XIII 
compare the results for the adaptive simulations of the 
models using the two different state vectors. Tables 
ieeeand XV show the tracking performance of the filter in 
the adaptive mode. As in Chapter III the tabulated data 
represent the time-averaged means of estimation error and 
the RMS value of the deviation of the means of estimation 
auper <£rom their time average. 
Wee SIMULATION RESULTS FOR THE CORRELATED RANDOM 

ACCELERATION MODEL SWITCH-ON-RANGE ADAPTIVE FILTER 

The correlated random acceleration model was also 
tested in an adaptive mode of operation against trajectories 
9 - 14. Table XVI shows this filter's position prediction 


performance results. Tables XVII - XIX show the filter's 
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Table XII 
SWITCH-ON-RANGE ADAPTIVE FILTER SIMULATION 


POR THE CONSDANT=-VELOGLLY MODEL 


TRAJECTORY SHIT SCARE 
9 90.58 5.55 

10 G0 o7 12.98 

11 71.58 11.39 

WP 63.95 Saas 
13 67.37 9.98 

14 74.67 9.65 
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Table XIII 
COMPARISON OF RESULTS FOR THE SPHERICAL AND 
CARTESIAN COORDINATE STATE VECTOR CONSTANT - 
VELOC Py SHOvEL owl TGH-ON-RANGE ADAPTIVE 


Pept ERS 

clare VEGTOR 
trade CTORY SPHERICAL CARTESIAN 
oHIT * OCARE 6HIT % 9CARE 
9 SIS St: So 1.38 A) 
10 607317 Arar: 47.20 10.60 
11 : fa es 2 1 ey 60.10 14.20 
eZ OS:29 5 ls.45 B20 13.40 
PS OF Si oe GU. 20 P5200 
14 74.67 OF0W, 54.70 oe) 
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tracking performance in the adaptive mode. As in the pre- 
ceding section, the tabulated data represent the time- 
averaged means of estimation error and the RMS value of the 


ewmration of the means of estimation error from their time 


average. 
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Table XVI 
SWITCH-ON-RANGE ADAPTIVE FILTER SIMULATION 


FOR THE CORRELATED RANDOM ACCELERATION MODEL 
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V. SUMMARY AND CONCLUSIONS 


A. SUMMARY OF TEST RESULTS 
1. Summary of Simulation Results for the Two Models 

The results presented in Chapter III show that 
against non-maneuvering target profiles the optimized 
constant-velocity models provided superior prediction 
performance. When the bandwidth of the filter is increased, 
adequate position prediction performance against maneuvering 
meee tS 215 Obtained Dut at the expense of performance 
degradation against the non-maneuvering profiles. The 
optimized correlated random acceleration model showed 
only a slight improvement against the maneuvering profiles 
considered, but produced results inferior to those obtained 
from the constant-velocity model against the non-maneuver- 
moemuare et profiles. The performance of the constant- 
velocity model against long-range maneuvering profiles 
was Significantly lower than for the short-range maneuvering 
profiles. This is to be expected since the long-range 
profiles require that the predictions be made over a longer 
period of time during which the target may maneuver 
unexpectedly and significantly. The tabulated data 
presenting the tracking performance of the filters 
developed in Chapter II, show that the position estima- 
tion performance for both filters is approximately the 


same. The correlated random acceleration model gave only 
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memeatiy better rate estimation performance than did the 
constant-velocity model. Against the target profiles 
Semsiaered the filter models derived from a spherical 
Soeordainate state vector provided significantly better 
overall predicted position results than the filter models 
derived from a Cartesian Coordinate state vector. 


2. Summary of Simulation Results for the Adaptive 
1 eens a 


Miceresmutserepontcde an Chapter LY for the 
Simulation of the switch-on-range adaptive filters indi- 
cated that the constant-velocity model derived from a 
femenrcal coordinate state vector provided significantly 
better predicted position performance than the constant- 
velocity model filter derived from a Cartesian Coordinate 
eete vector against all target profiles considered. 
Further, the results indicate that of the models derived 
meodme spherical coordinate state vector, the constant- 
velocity model gave the best overall predicted position 


pommermance . 


Prec ONCLUSIONS 

A filter derived from a spherical coordinate state 
vector with assumed linear system dynamics produces pre- 
dicted position results potentially superior to those 
obtained from a filter derived from a Cartesian Coordinate 
Pate vector. The spherical coordinate state vector 


constant-velocity model provides better overall predicted 
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position results than the correlated random acceleration 
model derived from a spherical coordinate state vector. 
The reason for this is perhaps due to the manner in which 
the predicted position, X(K+N/K) is generated. For both 
models considered, the predicted position x (K+N/K) 1s formed 
by raising the state transition matrix o(K,K-1), to be 
the Nee power as shown by Equation (C.2). In the case of 
the constant-velocity model such a method uses the current 
position and velocity estimates to generate the predicted 
position X (K+N/K). The correlated random acceleration 
model uses the current position, velocity and acceleration 
estimates to generate the predicted position. Thus, any 
error in the estimate of the target's acceleration state 
causes significant errors in the predicted position when 
propagated over the prediction interval. These accelera- 
tion induced prediction errors are not present when only 
the current estimates of position and velocity are used 
Bemeieethe Case of the constant-velocity model. 

A filter model derived from a spherical coordinate state 
vector can produce results potentially superior to those 
obtained from the corresponding Cartesian Coordinate state 


vector model in an adaptive filtering environment. 


C. SUGGESTIONS FOR FURTHER INVESTIGATION 
Since the constant-velocity model produced the best 
overall predicted position results, a complete analysis 


@t the sensitivity of this model's predicted position 
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performance to the statistics assumed for the measurement 
noise process v(K) and the sampling period T should be 
carried out. Additional sensitivity studies should 

also focus on allowing the diagonal elements of Q(K) to 
vary independently, thus removing the constraint that 

9 iS 

For the correlated random acceleration model, one area 
that warrants future investigation is the use of only the 
weetreron and velocity estimates of the filter in generating 
the predicted position, and comparing the results with those 
obtained from the constant-velocity model filter. This 
would reduce the errors in the predicted position due to 
emgors im the current estimate of the target's accelera- 
iewon State. 

Both models should be tested against a greater variety 
Gt target profiles, especially targets that present crossing 
mmeorales in which large bearing and elevation rates are 
generated, and targets that execute maneuvers while in 
itevel flight. 

Finally, more sophisticated adaptive filtering schemes 
should be developed and investigated for the constant- 
velocity model. A prime candidate for such an investigation 
is the residual testing method reported by Price and 


Brown [Ref. 11]. 


63 





APPENDIX A 


MONTE CARLO SIMULATION PROGRAM 


A.1 PROGRAM LISTING AND DESCRIPTION 

The Monte Carlo Simulation Program used in this thesis 
was written by Dr. D. E. Kirk of the Naval Postgraduate 
School. It 1s extremely flexible and well documented to 
Mesist the potential user in data preparation and input, 
as well as the various options available. It should be 
pointed out that SUBROUTINE KILL and all statements in the 
main program that deal with miss distance data were written 
Seeoeitically for this thesis. They may be deleted with no 
loss of program flexibility or functional operation. 

Computer system requirements vary with the simulation 
and the user's desires. All the computer work done for 
this thesis was performed on the IBM SYSTEM 360/67 at the 
feaval Postgraduate School. The program as currently 
dimensioned requires 520K bytes of storage. When redimen- 
fwomed for the 60 data point trajectories, 320K bytes of 
storage are required. A 100-member ensemble simulation of 
peeeiniear Sixth-order filter with no on-line computation 
meauired, the trajectory read in, miss distance information 
ered, and all printing and plotting options invoked 
requires approximately 1.5 minutes of CPU time for execution. 
An additional 30 seconds (approximately) are required for 


the compile/link process. For the work done for this thesis 
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the program was pre-compiled (FORTRAN IV Level G) and 
stored and executed from a data cell. The link process 

in this case required approximately two seconds. If the 
same linear, sixth-order filter is again simulated but the 
printed and plotted output is suppressed (retaining only 
the miss distance data) 59 seconds of CPU time is required 
to execute the simulation. In the case of the linear 
ninth-order filters that were simulated, the storage 
requirements remain unchanged, as do the time require- 
ments for the compile/link process. A 100-member ensemble 
simulation of the linear ninth-order filter with no on- 
line computation required, miss distance information 
desired, the trajectory read in and all printing and plotting 
options invoked requires approximately 3.0 minutes of CPU 
time for execution. If only the miss distance information 
tsedesired and all other output is suppressed, approxi- 
Weepeky 2,55 minutes of CPU time is Required, COVexccutLe 


the simulation. 
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A.2 LISTING OF KEY VARIABLES USED IN THE MONTE CARLO 
rei. <———_—<"  —"  °.. 


ROGRAM 


AMISS 


COVW(m,m) 


Een ni ) 
ERR(N,nsam) 


G(n,m) 
GAMMA (n,m) 


GAMMAS (n,m) 


GK (n,m) 
GKS (n,m,nsam) 


ecm, n ) 


HIT 


HS (m,n) 


iE R 
M 


MISSD(nens ,nsam) 


N 

NENS 
NSAM 

PHI (nn) 


average number of misses over the entire 
track. (scalar). 


covariance of the system random forcing 
input matrix (double precision). 


identity matrix (double precision). 


Sample mean of estimation error vector 
(sane le preersaen) . 


Sample gain matrix (double precision). 


System random forcing input distribution 
matrix (double precision). 


System random forcing input distribution 
matrix (single precision). 


Sample gain matrix (Single precision). 
filter gain schedule (single precision). 


system observation matrix (double 
precision). 


average number of hits over the entire 
track (scalar). 


System observation matrix (single 
Dieee ls ion. 


ensemble member index. 
number of measurements. 


sample miss distance matrix (double 
precision). 


system order. 
number of members in the ensemble. 
number of time points to be processed. 


State transition matrix (double 
DieeGis hon). 
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Peo (n,n) 


PKK(n,n) 


PKKMI1 (n,n) 


PKKM10 (n,n) 


Eiko (n,n,nsam) 


Q(m,m) 


R(m,m) 


SCARE 


SIGV(m) 


SI GW (m) 


SIGXZ (n) 


V (mn) 


VAR(n,n,nsam) 


W(m) 


X (n) 
XHATZ (n) 


XHKK (n) 


XHKKM1 (n) 


Sire meqansitione matrix (Single 
precision). 


theoretical covariance of estimation 
error matrix P(K/K) (double precision). 


one-step prediction error covariance 
matrix P(K+1/K) (double precision). 


initial value of the covariance of 
ecelmatiomrecinal Matrix GW 1) (double 
DRC eas LON). 


sample covariance of estimation error 
matrix P(K/K) (single precision). 


System random forcing input covariance 
matrix (double precision). 


Govdtlanecenor Measurement molse Matrix 
(double precision). 


average number of Scares over the entire 
rae ke 


standard deviation of measurement noise 
Meccor s(sinale reels ion) . 


Standard deviation of system random 
forcing input vector (Single precision). 


Standard deviation of system initial 
condition vector (single precision). 


measurement noise vector (single 
precision). 


Sculd ecOmaltance On eStiMation error 
matrix (single precision). 


system random forcing input vector 
(Single precision). 


state vector (single precision). 


PPE eomVvalUueTOn the state estimate 
vector =(C//= Iee(sanele precision). 


. State estimate vector x (K/K) (single 


Drees lon). 


one-step state prediction vector x (K/K- 1) 
(cine precision): 
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XM(n,nsam) 


XS (n,nsam) 


XZMEAN (n) 


mean state trajectory over the ensemble 
(Single precision). 


reference trajectory (single precision). 


mean of initial condition vector 
(Single precision). 
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A.3 IBM SYSTEM 360/67 SOURCE LIBRARY SUBROUTINES CALLED 


BY THE MONT IC M 
GAUSS3 matrix inversion by the method of GAUSS. 
OVFLOW utility program required by the random 
number generator. 
EOP generates two-dimensional line printer 
pOMLONE S 5 
SNORM random number generator that gives 


normally distributed, shuffled pseudo- 
random numbers. 
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A.4 MONTE CARLO SIMULATION PROGRAM FLOWCHART 


START 












READ 
o(K,K-1) 
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READ = 1 = 1 | READ 
P(K) Q(K) Pr! (K) Q(K) 
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7 0 
READ 
P(0/-1) 
WRITE 
IP (0/-1) 
# 0 
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READ 
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WRITE MEAN 
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COMPUTE 
SPAT oEl CS FOR 
THE MONTE 
CARLO 













WRITE G(K) 
P(K/K) FOR 
EACH K 







WRITE SAMPL 
MEANS AND 
VAR. OF 

Bye UUR eb On 










IEPAA) 










PLOT 
Pahe/ yas, . 
rs K 






PLOT MEAN 
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MEAN OF 
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Vike Or 
ESTIMATION 







Hee 





COMPUTE MISS 
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PLOT MEAN 
OF MISS 
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APPENDIX B 
DESt ee URACECTORIES 


Bet. DESCRIPTION OF TEST TRAJECTORY PROFILES 


a. Trajectory 1; short-range, low speed, non- 
maneuvering 


Trajectory 1 commences at a range of 3850 yards 
in a 30 degree dive with an initial velocity of -230 
yds/sec (the negative velocity indicates motion toward 
the origin). The trajectory is terminated at a final 
range of 457 yards. No maneuver is executed. Sixty 
data points were taken at a 4-Hz rate. 


b. Trajectory 2; short-range, low speed, moderate 
level maneuver 


Trajectory 2 commences at a range of 3850 yards 
in a 30 degree dive with an initial velocity of -230 
yds/sec. At 3450 yards the target executes a 2.0 g sin’ 
planar maneuver lasting 13.5 sec. The maneuver is completed 
at the final sample range of 465 yards. Sixty data points 
were taken at a 4-Hz rate. 


wer ibaTecctOn eo One-range, low Speed, non- 
maneuvering 


Trajectory 3 commences at a range of 9225 yards 
imma su ederree dive with an initial velocity of -230 
yds/sec. No maneuver is executed. The trajectory is 
terminated at a final range of 5832 yards. Sixty data 


points were taken at a 4-Hz rate. 


isi 





d. Trajectory 4; long-range, low speed, moderate 
level maneuver 


Trajectory 4 commences at a range of 9225 yards 
in a 30 degree dive with an initial velocity of -230 
Mde/sec. At 8822 yards, the target executes a 2.0 g pam 
mlanar maneuver lasting 13.5 sec. The maneuver is com- 
pleted at the final sample range of 5840 yards. Sixty 
data points were taken at a 4-Hz rate. 


e. Trajectory 5; medium-range, high speed, evasive 
level maneuver 


Trajectory 5 commences at a range of 4500 yards 
gna 20 degree dive with an initial velocity of -300 
yds/sec. At 3900 yards the target executes a 6.0 g enna 
iranar maneuver Lasting 10 sec. The maneuver is completed 
at a range of 914 yards and the trajectory is terminated 
at a final range of 90 yards. Sixty data points were 
taken at a 4-Hz rate. 


f. Trajectory 6; long-range, high speed, evasive 
evel maneuver 


Trajectory 6 commences at a range of 8000 yards 
ied) 20 aegree dive with an initial velocity of -300 
yds/sec. At 7500 yards the target executes a 6.0 g sin? 
planar maneuver lasting 10 sec. The maneuver is completed 
at a range of 4340 and the trajectory is terminated at a 
mange Of 5590 yards. Sixty data points were taken at a 


4-Hz rate. 


ro 





g. Trajectory 7; medium-range, high speed, non- 
maneuvering 


Trajectory 7 commences at a range of 5000 yards 
in a 20 degree dive with an initial velocity of -300 
maey,sec. Nommaneuyer 1s “executed. The trajectory is 
terminated at a final range of 575 yards. Sixty data 
peants were taken at a 4-Hz rate. 


h. Trajectory 8; long-range, high speed, non- 
Maneuvering 


Trajectory 8 commences at a range of 8000 yards 
imide 20 degree diye with an initial velocity of -300 
yds/sec. No maneuver is executed. The trajectory is 
terminated at a final range of 3575 yards. Sixty data 
points were taken at a 4-Hz rate. 


1. Trajectory 9; long-range, low speed, non- 
maneuvering 


DGayeetony commences at a range of 3000 yards 
in 20 degree dive with an initial velocity of -200 yds/sec. 
No maneuver is executed. The trajectory is terminated 
at the final sample range of 50 yards. One-hundred-sixty 
data points were taken at a 4-Hz rate. 


j. Trajectory 10; long-range, low speed, moderate 
level maneuver 


Trajectory 10 commences at a range of 8000 yards 
in a 20 degree dive with an initial velocity of -200 yds/sec. 
At 4000 yards the target executes a 2.0 g sin’ planar 
maneuver lasting 13.5 sec. The maneuver 1S completed at a 


range of 1200 yards and the trajectory is terminated at 


io 








the final sample range of 50 yards. One-hundred-sixty 
data points were taken at a 4-Hz rate. 


k. Trajectory 11; long-range, low speed, moderate 
evel maneuver 


Trajectory 11 commences at a range of 8000 yards 
in a 20 degree dive with an initial velocity of -200 yds/sec. 
mez 700 yards the target executes a 2.0 g sin? planar 
maneuver lasting 13.5 sec. The maneuver is completed at 
the final sample range of 50 yards. One-hundred-sixty 
data points were taken at a 4-Hz rate. 


1. Trajectory 12; long-range, low speed, moderate 
level maneuver 


Trajectory 12 commences at a range of 8000 yards 
in a 20 degree dive with an initial velocity of -200 yds/ 
sec. At 3333 yards the target executes a 2.0 g sin’ 
planar maneuver lasting 13.5 sec. The maneuver is completed 
at a range of 500 yards and the trajectory is terminated 
at the finai sample range of 50 yards. One-hundred-sixty 


data points were taken at a 4-Hz rate. 


m. Trajectory 13; long-range, low speed, moderate 
level maneuver 


Trajectory 13 commences at a range of 8000 yards 
in a 20 degree dive with an initial velocity of -200 yds/ 
sec. At 2700 yards the target executes a 2.5 g sin’ 
planar maneuver lasting 10 sec. The maneuver is completed 


at the final sample range of 50 yards. One-hundred-sixty 


data points were taken at a 4-Hz rate. 











n. Trajectory 14; long-range, low speed, moderate 


level maneuver ©... 

Trajectory 14 commences at a range of 8000 yards 
in a 20 degree dive with an initial velocity of -200 yds/ 
Bece, At 4000 yards the target executes a 2.5 g sin’ 
milanar Maneuver lasting 10 sec. The maneuver is completed 
meeaetange of 6000 yards and the trajectory is terminated 
at the final sample range of 50 yards. One-hundred-sixty 


data points were taken at a 4-Hz rate. 


B.2 TRAJECTORY GENERATING PROGRAM 

Giewtnaiecetories of Section B.1 are the products of a 
Pracking generating program written by Parr. A complete 
description of its use, development of the planar sin 
maneuver, program listing and several examples may be 


fiomme in Parr’s work jRef. 4]. 
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APPENDIX C 
Tibet eoter PREDICTION PROBLEM 


The ultimate test of any fire control system is its 
ability to destroy an engaged target. To this end, the 
determining factor is the accuracy with which it predicts 
the future position of the target, since launcher laying 
orders are derived from knowledge of the target's present 
Miceauture States. This of course, neglects the response 
characteristics of the launcher and the weapon (projectile) 
bagdastics, both of which are assumed here to have 
meeligible effects on the end result. 

The interval over which the predictions are made is 
derived £rom a straight line approximation of the plot of 
Paeempectile Time of Flight (TOF) vs. the target's present 
meme tange, £0r a 5/54 CAL. projectile. For the inter- 
wiein Question, {TOF: TOFe(0.0,10.0)} the straight line 


approximation takes on the value 
TOF = 1.43 x Ly (sec/yd) x (slant range) (C.1) 


The TOF interval, (0.0,10.0) was based on the assumption 
jiac the maximum effective range of the 5/54 CAL gun 
against an air target is 7000 yds, corresponding to a 


TOF of approximately 10.0 seconds. 
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In the absence of any deterministic inputs to the 
Eaten yetne tareet's future poSition predicted over N 


samples is given by 
X(K#N/K) = g°X(K/K) (C. 2) 


Raising the state transition matrix > (K,K+1) to a decimal 
power PEcScieo SOMe meOmMpUGcatlonaledsifttenlties, so N is 
formed by computing the TOF from X, (K/K) (the current 
estimate of the target's slant range) then truncating this 
value to give the integer N. For large N the difficulty 
Sreraising a (9 x 9) matrix to the ne® power SeimeexdS ts. 
Meee, however, that the assumption of linear plant dynamics 
decouples the filter in each of three coordinate directions. 


Partitioning and decoupling we have for the systen, 


ice 0 x (K)] PT, 0 0 | Jw, CK) 
X(K+1) = [0 Qo. XoCK)P+)0 Fo 98 | jw, Gd) (Cc. 3) 
| 
0 0 053 lS E oat es 
where 
x, (K) = [r,r,r]" w(K) = [w, (K)] 
x5(K) = [6,6,6]° w(K) = [w,(K)] (C.4) 
x3(K) = [0,6,6]° = w(K) = [w,(K)] 


IS) i 





Further, note that for each of the models discussed in 
this thesis, each of the Oa3 iis are tower tri- 
angular matrices. Following the method given by Parr 
ier, 4], the elements of the bs 1 = 1,2,3 are operated 
on individually and the resulting by 1 = 1,2,3 are 
Pyoined to form the system state transition matrix 

o. (K#N,K). The intermediate result follows quickly and is 


Shown on the following page. 
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Finally, adjoining the ee gives 


: 0 
oN (KN, K) = 65, 0 (C.6) 
N 
: 35 





The result shown in Equation (C.5) is easily programmed 
and involves simple manipulations without the need for 
extensive matrix multiplication. Note that the result 
Shown in Equation (C.5) 1s Baaett POU(EOr theysinele 


dimension) 
r(K#N) = r(K) + £(K)[(N-K)T] + 3 TCK)[(N-K)T]® (C7) 


For the case of the constant-velocity model, o,,(1,3), 

o.;, (2,3), and $,,(3,j) ieee. SG ame Set £o. zero. The 

Piecatction ts then easily found by using Equation (C.2). 
The magnitude of the N-step prediction error is 

defined to be the magnitude of a position vector, say 

——_—* « ° 

ERRNET, drawn between the eStimate at time t = (K+N)T, 


mukeN/K), and the reference trajectory, x(K+N) 


| |ERRNET|| = ||x(K+N) - X(K+N/K) || (OLD 


Note that no reference is made to the direction of the 
error vector since the previously defined kill region is a 
sphere of constant radius centered at the reference 


position x(K+N). 
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iiecwpminin simeoteseeides insehe following three 


regions 
1) | |ERRNET | | < 8.0 yds nai ze 
ee 
Zee Om ERRNED eer. 0 yds. scare 
— See 
3) | |ERRNET|] > 16.0 yds miss 


To avoid the necessity for coordinate conversion, the 
error vector is determined by a simple rotation of the 
polar plane followed by applications of the law of coSines. 
Consider the situation depicted in Figure 4, where points 
a and b represent X (K+N) and x(K+N/K) (position only). 


Then from Equation (C.8) 


||ERRNET|| = |[ab]] = | [bal | (C.9) 


If the polar plane in Figure 4 is rotated through some 
angle 6 about an axis constructed by drawing a line through 
the origin connecting the points (090,270) (keeping Ty 


ama r. fixed) Figure 5 is obtained. The angle of rotation 


2 


isedetined as 


Bs min(®, ,9,) (C.10) 


Point now lies on the rotated polar plane, the position 
of point b remains fixed, and point c represents the fixed 
range projection of point b onto the polar plane. Some 


rather obvious geometric relationships follow immediately; 
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A® = /boc = |8, - 8 


AR = jad| = |r, - r,| (G1) 
=i 7 
lod] = max (ry 51.) 

The problem then remains of finding |ab]. Figure 6 shows 


the projection onto the rotated polar plane. Note that 


triangle cod is an isosceles triangle, then 


_ 180 - Ad (Ca) 


Applying the law of cosines lac | is given by 


Jac] = {(lad])? + [(2.0) (sing®,) (Jod|) 17 - 





) (6413) 
(2.0) (fad]) (2.0) (sing®,) (fod]) cos/cdo) }1/? 
[bc | iS given by 
[be | = (2.0) ([od|) (sind?) (Clr) 


Again applying the law of cosines and noting that 


Pa | 
/bced = /bco, |ab/| is given by 


eZ 


fab] = {|bc|? + Jac|? - (2.0)(|be]) Clac])cos (/bed)} Gos 


Thus, the magnitude of the N-step prediction error has been 
determined without the need for any coordinate conversion. 


The validity and accuracy of this means of determining the 
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error magnitude was verified by conversion of x(K+N) and 
and x (K#N/K) to the Cartesian Coordinate system and deter- 
Mining the magnitude of the error vector by standard means. 
In every case investigated, the results given by the two 
different methods were in complete agreement. Again note 
that no attempt is made to determine the direction of the 
Seon VECTOT rab PiMecmEnc hectons Of Interest Speci iy 
concentric spheres of constant radii centered at X(KtN). 
The directionality of the error vector is eliminated in 
SUBROUTINE KILL through the use of the absolute values of 
the differences in bearing and elevation, and by using the 


maximum value of the actual and estimated ranges. 
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